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■ Spin-orbit corrections to superexchange are calculated using the method of 

Moriya [T. Moriya, Phys. Rev. 120 91, (I960)] for two of the insulating parent 
compounds of the cuprate superconductors: (1) La2- x Nd x Cu04 where the CuC>6 
octahedra forming each Cu-0 layer are tilted in staggered fashion about an axis 
£S) ' which depends on x and temperature; and (2) YBa2Cu306+x (x < 0.4) where the 

Cu-0 layers form CUO2-Y-CUO2 bilayers in which the in-plane 2 ~ ions are dis- 
placed uniformly towards the Y 3+ layer. For (1) a simple formula is derived for the 
weak ferromagnetic moment in each Cu-0 layer as a function of the tilting axis and 
magnitude. For (2) it is shown that the anisotropic corrections to superexchange are 
different from what has previously been assumed. For the correct spin Hamiltonian 
a classical Neel state in which the Cu spins are lying in the plane is unstable in a 
single Cu-0 layer, but when a bilayer is considered there is a critical value of the 
interlayer exchange coupling which stabilizes this state. For both cases (1) and (2) 
spin-wave spectra are calculated and shown to compare favorably with experiment. 
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I. INTRODUCTION 



Spin-orbit (SO) coupling causes electron spins to precess as they move through the elec- 
tric field of a crystal lattice. Within the tight-binding approximation this precession appears 
as a small spin rotation which occurs whenever an electron tunnels between two Wannier 
orbitals. As first shown by Moriya |[| this rotation can have important consequences in 
antiferromagnetic (AFM) Mott insulators; when it is included in Anderson's calculation 
of superexchange then anisotropic corrections to the otherwise isotropic effective spin Hamil- 
tonian are generated. These corrections, known as Dzyaloshinski- Moriya (DM) interactions 
lift any ground state degeneracy associated with rotational invariance in spin space 
and are responsible for such effects as weak ferromagnetism and the existence of spin-wave 
anisotropy gaps. 

The subject of this paper is the DM interactions which exist in the distorted Cu-0 
layers of the insulating AFM parent phases of the cuprate superconductors. We will be con- 
cerned both with the microscopic origin of these interactions and their physical consequences. 
Probably the best known example of the effects of DM interactions in the cuprates occurs in 
La2Cu04 ||. In this material a structural phase transition occurs from a high temperature 
tetragonal (HTT) phase (space group 14/mmm) to a low temperature orthorhombic (LTO) 
phase (space group Bmab) . In the LTO phase the Cu06 octahedra forming each Cu-0 layer 
tilt in a staggered pattern about the (110) axis and this distortion results in DM interactions 
which induce a weak ferromagnetic moment in each layer M. Another example of the effects 
of DM interactions in the cuprates is the easy-plane anisotropy which has been observed in 
the spin- wave spectrum of the AFM insulating phase of YBa 2 Cu30 6+x (x < 0.4) |J. In this 
case the relevant structural feature is that the Cu-0 layers form CUO2-Y-CUO2 bilayers in 
which the negatively charged in-plane O 2- ions are uniformly buckled towards the positively 
charged Y 3+ layer. 

One motivation for the present work is that in a recent paper a 'one-band' description 
of SO coupling in both insulating and doped La 2 Cu0 4 in the presence of various tilting 
distortions was studied using the Hamiltonian 

#flux = -t £ {e^ a clc ]a + h.c) + C/$> iT n a . (1.1) 

(ij) * 
a 

Here c\ a creates an electron with spin a at site i, n ia is the corresponding number operator, 
~ (— 1)^+*) 0.050, (fii^+y = —<fii,i+x where 9 is the octahedral tilt angle, and in the 
exponent a = +(— )l/2 for up (down) spins. The x and y components of site i are denoted Xi 
and yi, and x and y are unit vectors. Hamiltonian Ql.lD describes a correlated tight-binding 
band of electrons which move in a background of staggered flux where up and down spin 
electrons have opposite charge. This flux is simply the Berry's phase associated with spin 
precession about the z-axis in spin space. In W\ the effective Hamiltonian describing the large 



U/t limit of (1.1) was derived and at half-filling the classical spin ground state was found to 
have no weak ferromagnetic moment. It is shown here that this result is not in conflict with 
the experimental observation of such a moment in the LTO phase of La2Cu04. Hamiltonian 
( |1.1|) is a correct description of electrons in a tilted Cu-0 layer when the system is viewed 
in the appropriate site- dependent coordinate system in spin space. When local rotations are 



performed to transform the system back to the physical spin-space coordinate system a weak 
ferromagnetic moment which agrees with experiment appears. 

A similar approach to weak ferromagnetism has been discussed recently by Shekhtman, 
Entin-Wohlman and Aharony (SEA) || who also studied the DM interactions induced by 
tilting distortions in L^CuO^ SEA were able to show that, quite generally, the DM in- 
teractions present on a single Cu-O-Cu bond are isotropic when the bond is viewed in the 
appropriate local coordinate system in spin space. This observation led them to the interest- 
ing conclusion that any physical anisotropy (i.e., anisotropy which cannot be 'gauged away' 
by local rotations) must arise from the frustration of these bonds. By applying this idea to 
the LTO phase of La2Cu04 SEA were able to successfully account for the observed weak 
ferromagnetism in this material. In this paper we calculate the DM interactions which occur 
in the presence of tilting distortions in which the CuC>6 octahedra can tilt about any axis, 
not just (110). Such general tilting distortions may have physical relevance because they 
describe at least the average structure of La2_ x Nd 1 .Cu04 when x ^ 0.5 P,|TU|]. The main 



result of our analysis is that regardless of the octahedral tilt axis the effective Hamiltonian 
describing a single layer can always be transformed into ( |1.1[ ), to lowest order in 9. However, 
the coordinate system in spin space in which ( [Lip holds changes as the tilt axis changes so 
that the ratio of the weak ferromagnetic moment to the size of the tilting distortion depends 
on the tilt axis in a simple way which we derive in Section IV. These results provide a po- 
tential experimental test of the 'anisotropy through frustration' idea of SEA [§] and justifies 
the use of (TO) in M. 



A second motivation for this work is the puzzling observation by Coffey, Rice and Zhang 
TTJ that the DM interactions in a single buckled Cu-0 layer of YBa 2 Cu306+ x tend to 
stabilize an incommensurate spiral spin configuration, while neutron scattering experiments 
see no sign of such a spiral ||. It is shown here that if an effective spin Hamiltonian for a 
single Cu-0 layer in YBa 2 Cu30 6+2: is derived using the same method as for La2_ x Nd. r Cu04 
then a classical Neel state with spins lying in the xy plane is indeed unstable; complex 
frequencies corresponding to exponentially growing unstable modes appear in the classical 
linearized spin-wave spectrum. However, when instead of a single Cu-0 layer a Cu02-Y- 
Cu0 2 bilayer is considered these complex frequencies disappear for a critical value of the 
interlayer coupling Jf 2 — 1 x 10~ 3 J, well below the lower limit on J 12 set by experiment ||, 
and the Neel state becomes stable. This stabilization occurs because the DM interactions in 
the upper and lower planes favor spirals with opposite senses (a consequence of the inversion 
symmetry of the YBa 2 Cu306+a: unit cell) and this spiraling is frustrated by the interlayer 
coupling. We also find that when J i2 is large enough not only does the Neel state become 
stable, but the spin-wave spectrum shows an in-plane gapless mode and an out-of-plane 
gapped mode in agreement with experiment 0. This result shows that the 'easy-plane' 
anisotropy in YBa 2 Cu306 +x in fact arises from frustration of DM interactions in accordance 
with the general principle of SEA ||. 

This paper is organized as follows. In Section II the structural and magnetic properties 
of La 2 _ x Nda;Cu04 and YBa2Cus06+x which are relevant for this paper are discussed. The 
SO modification of superexchange due to the structural distortions in these materials is 
calculated in Section III and the resulting classical ground state and spin-wave excitation 
spectra are presented in Section IV and V, respectively. Finally, Section VI summarizes the 
conclusions of the paper. 



II. RELEVANT EXPERIMENTAL FACTS 



A. La 2 _ a; Nd r Cu04 

The La 2 _ a: Nd x Cu04 system shows a rich structural phase diagram as a function of x and 
temperature. As mentioned above when x = the material undergoes a phase transition 
from the HTT phase to the LTO phase. It is also known that when enough Nd is doped into 
the system (x > 0.5) the material crystallizes into the T' structure of pure Nd 2 Cu04 [|12 



However, for smaller Nd concentration (x < 0.4) there is still an HTT — > LTO transition, 
and as the temperature is lowered further there is a second transition into a structural 
phase with space group Pccn PJT0|] . In both the LTO and Pccn phases the Cu0 6 octahedra 



forming each Cu-0 layer tilt in a staggered fashion through an angle #(~ 0.1) about first the 
(cosx, sinx, 0) and then the (sinx, cosx, 0) axis in successive Cu-0 layers where % — 7r /4 in 
the LTO phase and 0<x<7r/4in the Pccn phase. The case \ — corresponds to the low 
temperature tetragonal (LTT) phase (space group P4 2 jncrri) which occurs, for example, in 
the doped material Lai.88Bao.i2Cu04 fl~3|j . 

Neutron scattering measurements of the spin structure factor of insulating La 2 Cu0 4 
and subsequent theoretical analysis have shown fairly conclusively that the spin degrees of 
freedom in this material are well described by an AFM Heisenberg model with exchange 



coupling J ~ 130 meV ||I4[. Although this model is adequate for describing most properties 
of La 2 Cu04 slight deviations from perfect isotropy have been observed experimentally. In 
particular: (i) Thio et al. and Cheong al. observed a first order spin-flop transition as a 
function of applied magnetic field perpendicular to the Cu-0 planes - this was interpreted 
as being due to the existence of weak ferromagnetic moments in each layer; and (ii) both 



neutron scattering |16| and AFM resonance measurements |17j have shown that the zone- 
center spin- waves in La 2 Cu04 are gapped with an in-plane gap of ~ 1.0 meV and an out-of- 
plane gap of ~ 2.5 meV. It was immediately realized by the groups which performed these 
measurements that these effects were manifestations of DM interactions. These interactions 
have since been calculated microscopically |TI]||[TS[]. The resulting spin Hamiltonian has an 
Ising-like anisotropy, which is responsible for the zone-center spin-wave gaps, and a ground 
state in which the spins lie nearly along the orthorhombic c axis except for a slight cant out 
of the Cu-0 plane which gives each layer a weak ferromagnetic moment. The canting angle 
is roughly Q w f ~ 0.005 and so the weak ferromagnetic moment is ~ 0.003/i# per Cu site. 



B. YBa 2 Cu 3 6+a; 

For x < 0.4 the YBa 2 Cu306+a; system is tetragonal and AFM with the added oxygens 
going into the chains presumably at random. When x > 0.4 the oxygen ions partially 
order and this leads to a tetragonal — > orthorhombic transition. At the same doping the 
planes become metallic and the material becomes a superconductor. Here we are primarily 
concerned with the tetragonal insulating AFM phase. The key structural feature in this 
material is that the Cu-0 layers are not equally spaced as they are in the La system but 
instead form Cu0 2 -Y-Cu0 2 bilayers. The planes forming these bilayers are buckled with 
the in-plane negatively charged 2 ~ ions displaced uniformly towards the positively charged 



and is roughly 0.22 A |19| and so the CuO bonds make an angle of roughly 6 ~ 0.1 with 



the Cu-0 plane [TJ5|]. Thus the magnitude of the distortion is approximately the same as in 
La 2 Cu04. 

Previously the magnetic structure of YP^CuaOe+z has been modeled assuming each 



bilayer can be described by a Hamiltonian of the form [20 



H — X! X! {j*Sa,iSa,j + JxyiSa^aj + ^a,i^a,j)} + ^12 X! ' ^2,4 (2.1) 
a=l,2 ij i 

The spin-wave spectrum of ( |2.1|) has four branches: a gapless in-plane mode; a gapped out- 



of-plane mode with a gap of 4:S^j2(J xy — J) J; and two high-energy branches with gaps of 

4v / 25'V 'J Ji2- Neutron scattering has shown the existence of a gapless in-plane mode with 
spin-wave velocity ~1.0 eV-A (J ~ 150 meV) and a gapped out-of-plane mode with gap ~ 4 
meV H and so the low-energy spectrum of ( |2.1|) agrees with the experiment. At the same 
time the high-energy modes have not been observed for energy transfers up to 50 meV putting 
a lower limit on the interlayer coupling of J 12 ^ 0.0 IJ H^l]]. The microscopic basis of (|2.1|) 



is not as firm as that of the corresponding spin Hamiltonian for La2Cu04 and in Section 
V we will show that, in fact, when the same methods which have been used successfully to 
describe La2Cu04 are applied to YP^CusOe+z the resulting spin Hamiltonian is different. 
Nonetheless, it is possible to reproduce the experimentally observed spin-wave spectrum 
using the new Hamiltonian. 



III. SPIN-ORBIT CORRECTIONS TO 
SUPEREXCHANGE 

Superexchange occurs when two magnetic ions interact through their mutual overlap 
with an intermediate diamagnetic ion |J. Recently SEA || have used Moriya's method 
]T|] to derive a fairly general expression for the anisotropic corrections to superexchange 
due to SO coupling on a single such bond for spin-1/2. Their result, which corrects some 
omissions in an earlier calculation |llj as well as a numerical error in Moriya's original paper, 



has the interesting property that it can be related to an isotropic interaction by a unitary 
transformation ||. A similar expression for the anisotropic superexchange on a single bond 
was implicit in and in this section its derivation is sketched to show that in fact the 
expressions in J7| and || are the same and to establish notation for the rest of the paper. 
In the presence of SO coupling a single Cu-O-Cu bond is described by the Hamiltonian 



-ffcuocu = ^2{d\ a (t pd 5 a p + i\i ■ a a/ 3)Pp + 4a(WW + ^2 • Oap)d(j + /i.e.} 

a/3 

+U d d(n <m n dll + n m n d2[ ) + A dp Y,plp a (3.1) 



The vacuum state for ( |3.1| ) is one in which the Cu 3d and O 2p shells are full. The operators 
d ia and then create holes with spin a at Cu site i and the oxygen a orbital, respectively. 
The (hole) energy splitting between Cu and O sites is A pd , the on-site correlation on the Cu 



site is Udd, and t p d, Ai and A2 are hopping integrals, the latter two arising from SO coupling 

The problem of deriving the effective spin Hamiltonian for (|3.1| ) simplifies upon applying 
the unitary transformation 



d 



/t 

1/3 



u 2/3 



E 

a 

E 



exp i tan 



exp i tan 




(3.2) 



a/3 



When is expressed in terms of the primed operators the result is 

#CuOCu = E {tpdAaPa + W^LPa + h ' C ) + U ( n ^ n U + ™2T n 2|) + A E^a ( 3 ' 3 ) 
a a 

where t p d = (tp d + A 2 ) 1 / 2 , (we will consider only the case where |Ai| = | A2 1 = A). The 
transformation (|3.2|) absorbs the spin precession induced by SO coupling into a redefinition 
of the local coordinate system in spin space. Such a transformation is possible because the 
bond is essentially one dimensional; i.e., there are no closed loops around which an electron 
can hop and acquire a finite spin precession which cannot be transformed away. 

The effective spin Hamiltonian resulting from ( |3.1| ) in the limit t p d ^> Udd, A p d can now 
be found using standard methods [p2]] . The result is if Bond = JS[ ■ S' 2 , with 

1 



J 



^ L dp 



A 



dp 



u dd 



(3.4) 



and S[ = (l/2)d'\ a a a i3d' 1/3 + 0(t/[A p d, U]). When the unitary transformation ( |3.2| ) is undone 
the final result is 



#Bond = J(SIS Z 2 + cos0(S 1 1, S 2 x + SfS y 2 ) + sin (j){S^S y 2 - S\Sl 



(3.5) 



where ~ |Ai — \2\ltpd is the angle through which an electron spin precesses when it hops 
from site 1 through the intermediate orbital to site 2, and where the z-axis in spin space has 
been chosen to be parallel to the precession axis Ai — A2. 

After some algebra it is possible to show that ( |3.5| ) is equivalent to the result obtained 
by SEA [§]. Hamiltonian ( |3.5|) was also derived in precisely the form given above in |F] but 
starting from a 'one-band' description in which the oxygen ions were not included explicitly. 
The equivalence between the one-band and three-band pictures at half filling is easy to 
understand. If the Hamiltonian for a single bond is given by 



#l-band = { C la(- t( W + *^12 * + h - C ] + U ( n l^ n H + ™2T™2|) 



(3.6) 



and if for a given t we choose U so that 4(t 2 + A 2 2 )/t/ = J and A i2 — t/t p d{\\ — A2), then 
similar arguments to those given above vield fl3.5l) when U t. In what follows we will 



adopt this 'one-band' approach and describe a given Cu-0 layer with a Hamiltonian of the 
form 



H = { C L (-tSocp + i\j ■ dap) Cp + h.c.} + Uj2 n it n n- ( 3 - 7 ) 

(ij) * 
a/3 

At half filling there exists an entire class of models with different t, U, and Xij values which 
yield the same effective spin Hamiltonian in the large U/t limit. Away from half filling this 
is no longer the case and it is necessary to perform a mapping from the full three-band 
model to an effective one-band model in order to find the appropriate parameters |7J. In 
what follows we are only concerned with the half-filled case and so for simplicity we take 
t — t p d^ 1.3 eV 



It remains to compute the Ay vectors. At half-filling each Cu ion in a Cu-0 layer has 
one hole in its 3c? shell which, in the absence of SO coupling, occupies the d x 2_ y 2 orbital. If 
on each Cu site a SO interaction 

H so =f3Y J L l -S l (3.8) 

i 

is included where Li and Si are the orbital angular momentum and the hole spin at site i, 
respectively, and (3 ~ 0.1 eV for Cu p^| , then higher crystal field levels - d xy , d xz , and d yz 



- are mixed into the lowest lying d x 2_ y 2 state. These admixtures then modify the hopping 
integrals and give rise to the A terms in ( |3.1|) . 

Moriya [0 derived an expression in second-order perturbation theory for the Xij vectors 
in (|3.7|) which for our purposes reads: 



? V r i/3 ^r(m,i\Li\Q,i) (m,j\Lj\0,j)* x 
Xij ~Xi-Xj~— *y(m, a) -± tji{m, a) j (3.9) 



where |m, i) is a crystal field split level in the absence of SO coupling labeled by m at site i 
with m = corresponding to d x 2_ y 2. The energy of the mth level is e m , and t i j{m 1 a) is the 
hopping matrix element in the absence of SO coupling between the Cu orbital m at site i 
and the oxygen a orbital between sites % and j. The relevant matrix elements of L are 

(x 2 - y 2 ,j\Lj\xz,j) = iy, 

(x 2 - y 2 ,j\Lj\yz,j) = -ix, (3.10) 

and so to evaluate ( p.9| ) only the matrix elements tij(m,cr) need to be determined. This is 
done below for La2- x Nd x .Cu04 and YBa2Cu306+ x - 



A. La 2 _ 1 Nd :E Cu04 

In the LTO and Pccn phases of La2_ x Nd :r Cu04 the average structure is one in which the 
CuOe octahedra forming each Cu-0 layer are rotated through an ang le (_l)(**+w)0 about 
the axis (cos x, sin x, 0). In general there are then two CuO bond angles: ±6sinx and 



Cu x z orbitals and oxygen a orbital in a typical bond in I^-^Nd^CuC^. For this bond the 
hopping from the rotated xz orbitals to the a orbital is, for small 9, 



ti,i+a(yz, o) - (3.11) 



and similarly 



t i)i+ y(xz,a) ~ 0, 



t i>i+y (yz,a) ~ (-l) Xi+yi V pd Jcos X . (3.12) 



where V p dn is the hopping amplitude between a Cu d xy orbital and an oxygen ir y orbital. 
The resulting A vectors are then 

A M+ x ~ (-l) Xl+w (Aicosx,A 2 sinx,0), 

A w ~ -(-ir +K (A 2 cosx,A 1 sin X ,0) (3.13) 

and 



Ai ~ 



Vpdirft 



A 2 ~ — 9, (3.14) 

A non-zero Ai, which arises to leading order from direct hopping between neighboring d 
orbitals, has been included in (|3.13|) for completeness and also to show that such a term 
does not affect the results that follow. If we take the reasonable values of (3 ~ 0.1 eV, 
e xz — e x 2_ y 2 ~ l.OeV and V p dn — 1.0 eV we obtain 

Ai ~ 

A 2 ~ 100 meV 9 (3.15) 



B. YBa 2 Cu 3 6+I 

The oxygen ion displacement in a buckled Cu0 2 -Y-Cu0 2 bilayer, like the tilting in 
La 2 _ 2 .Nd a; Cu04, results in a non-zero tij(xz,a) and tij(yz,a). Because the geometry is 
different (the d orbitals are not rotated, see Fig. 0b) the proportionality to 9 of these hopping 
integrals is larger than in La 2 _ x Nd :r Cu04 by a factor which we estimate to be roughly 4.8 



25], and so in this case 



t i4+ x:(xz, a) ~ 4.8 V pdn 9, 

t iii+ x(yz,a) ~ (3.16) 



and 

t iti+ y(xz,a) ~ 0, 



and the resulting A vectors are 



XgU=(0,A,0), 

XS^ = -(A,0,0) (3.18) 

where 

a ~ e (3.i9) 

(We will refer to the upper and lower layers as 1 and 2, respectively.) If the crystal field 
splitting in YBa 2 Cu306+:r and La^^NdzCuC^ are assumed to be the same then A123 — 5A214. 
An important point for what follows is that because of the inversion symmetry of the 

— * 

unit cell of YBa 2 Cu 3 6+a; the Ajj vectors in the lower layer are the opposite of those in the 
upper layer and so 

\ (2 ix = -(0,A,0), 

Xg ) + y = (A,0,0). (3.20) 



IV. MAGNETIC ANISOTROPY IN La^Nd^CuC^ 

In order to derive the effective spin Hamiltonian for La 2 _ a; Nd :r Cu04 it is necessary to 
consider the large U/t limit of (|3.7[) with A^- vectors given by ( |3.13|) when there is one 



electron per site. Before doing this it is useful to analyze the structure of this model by 
considering the motion it describes for a single electron. For this purpose it is natural to 
decompose Ay into 'frustrated' and 'unfrustrated' components as follows 

A w = (-l)* i+w {ai(cosx,sinx,0) + a 2 (cosx, - sinx, 0)} (4.1) 
Vt-y = -(-l) Xl+s/l {ai(cosx,sinx,0) -a 2 (cosx,-sinx,0)} (4.2) 

where a\ = (Ai + A 2 )/2 and a 2 = (Ai — A 2 )/2. To see why this decomposition is useful 
consider the hopping of an electron around a single plaquette according to ( |3.7| ) for (i) 
«i 7^ 0, a 2 = 0, and (ii) a% = 0, a 2 7^ 0. For case (i) (unfrustrated) the sign of the spin 
precession oscillates and no net precession occurs (Fig. |2|a) while for case (ii) (frustrated) the 
sign of the precession does not oscillate and there is a net precession of ~ 4a%/t (Fig. ||b). 
As emphasized by SEA || case (i) is special: When a 2 = it is possible to perform a 
unitary transformation which maps ( |3.7|) onto a precession-free (isotropic) Hamiltonian. 
The non-trivial physics is thus due to the frustrated precession and it is precisely this which 
is described by ( |1.1| ). 

The unitary transformation which eliminates the unfrustrated precession in ( |3.7| ) cor- 
responds to a local rotation in spin-space. For 0,a 2 /t 1 this unitary transformation is 
represented by the operator (see Appendix) 

U = exp [i y](-l)(*i+%)^(cos X , -sinx,0) ■ s\ (4.3) 



under which (|3.7| ) becomes, to leading order in a/t, 



UHU ] ~ £ [c\ a (-tS aP + tKj ■ a aP ) c jf3 + h.c.} + (4.4) 

(id) » 
a/3 

where A^ +i = (— l) Xi+?/,; ai(cos x, sin x, 0) and A^ + ^ = — A^ i+& . If a global rotation in spin 
space is then performed to bring the z-axis parallel to (cos x, sin x, 0) the result is ( |1.1| ) with 
0M+x = (-l) (a:i+2/i) tan- 1 (a; 1 /t) ~ (-1)^+^0.050 and <f) i>i+f = It follows that to 

leading order in 9 the energy spectrum of ( |3.7|) is independent of the tilt axis angle x an d 
depends only on the octahedral tilt angle 9. 

At half-filling and when U ^> t the effective spin Hamiltonian arising from ( |1 . 1|) is (up 
to an irrelevant constant) [|7] 

Hse = Jj2{ S tS- + coaMSfS? + SfS]) + smMSfS* - SfS*)}. (4.5) 

m 



The sin0,jjz • (Sj x Sj) term in (|4.5|) is minimized by a four sublattice state which only 
becomes stable when the magnitude of this term is larger than J [IS|. This is never the case 



here and so this term is completely frustrated, there is no spin canting, and the remaining 
cos 4>ij easy-axis lines the spins up parallel to the z direction. 

The classical ground state of ( f4.5|) has no weak ferromagnetic moment and so seems 
to disagree with experiment. However, this is no longer true in the physical spin-space 
coordinate system. Once the local rotations used to transform (|3.7|) into (14^ and also 



the global rotation which brought the z-axis in spin space parallel to (cosx, sinx, 0) are 
undone the spins which once were parallel to the z axis become nearly aligned along the 
(cosx, — sinx,0) direction except for a slight cant out of the Cu-0 plane with the canting 
angle given by 

Q wf = a2 5111 2 * „ 0.050 sin 2 X . (4.6) 

The spin configuration in the physical basis for a specific x between and 7r/4 is shown in 
Fig. |3], together with the frustrated and unfrustrated spin precession axes. In the LTO phase 
(x = 7r/4) (|4.6j ) agrees with experiment: the spins point nearly along the (110) direction (i.e., 
parallel to the orthorhombic c axis) and cant out of the plane through the angle 0„,/ ~ 0.005 
0]. As x decreases the ratio of the weak ferromagnetic moment to 9 also decreases and this 
decrease should be experimentally observable in the Pccn phase of La2- x Nd x Cu04. 

Given the classical ground state of ( |4.5| ) the linearized spin-wave spectrum can be cal- 
culated using standard methods [p6fl . The result is a twofold degenerate spectrum defined 



within the AFM zone with dispersion (here q is in units of one over the lattice spacing) 

uj(q) = 2JS ^4 — cos 2 <f)(cos q x + cos q y ) 2 — sin 2 <p(cos q x — cos q y ) 2 ^j (4.7) 

The two zone-center modes, which in the physical basis correspond to in-plane and out- 
of-plane excitations, have a gap of 4JS'sin0. For the parameters used above this is ~ 1.5 



meV and ~ 2.5 meV in La2Cu04 [|T^,|1^]. Note, however, that the theory presented here 
unambiguously predicts that the in-plane and out-of-plane spin-wave gaps should be equal, 
while experimentally these gaps have different values. It is probable that this discrepancy 
is due to dipolar interaction between spins, another source of anisotropic spin interactions 
which has not been included in this calculation. 



V. MAGNETIC ANISOTROPY IN YBa 2 Cu 3 6+I 

Before studying the properties of a Cu0 2 -Y-Cu0 2 bilayer it is useful to first consider 
the simple toy model shown in Fig. |j: Two coupled chains in the presence of a buckling 
distortion. This model has the advantage of simplicity and the basic physics is similar to 
that of the more complex two-dimensional bilayer. The spin Hamiltonian for this system 
can be written 

H (ID) = H (1D) + H (1D) + H (XD) (51) 

where 

H[ 1D) = ./X{^,ST.,,, + cos0(^ m + S^ )i+1 ) + sin^AV - 

(5.2) 

^2 ^ = J X/{ ^2,1^2,1+1 + COS 4>(^2,i^2,i+l + &2,i ^2,i+l ) — ^(S^iS^i+l ~ ^fjiS^i+l) } 

(5.3) 

and where we also include an isotropic interchain coupling 

H& D) = Ji 2 E 5 M-*V (5.4) 

i 

First consider the case J12 = 0. If we define the unitary operators 

UM = expii<f> y £x i S^. (5.5) 

then it is possible to transform away the DM interactions in the two chains as follows: 

U^fiH^Uttt) = J XX* • S w (5.6) 

i 

U 2 {-4>)Hl 1D) Ul{-(j>) = ■ S 2 ,i+i- (5-7) 

i 

The eigenvalues of Hi and H 2 are therefore the same as those of two uncoupled isotropic 
Heisenberg models. The corresponding eigenvectors are, however, different. For example, 
consider the classical ground state manifold of Hi. This manifold is infinitely degenerate 
because ( |5.6D , which is a unitary equivalent of Hi, is isotropic in spin space; but the ground 
states of Hi are not the same as the ground states of (|5.6|). Two extreme cases are (i) the 



to the x-axis and form a Neel state. For both cases the energy gain per site over having the 
spins lying uncanted in the yz plane is JS 2 sin and so these states are degenerate. 

The system changes qualitatively when the slightest interchain coupling is introduced. 
Figure ^ shows one of the degenerate classical ground states of the system when J12 = in 
which the spins spiral in the yz-plane. Note that the senses of the spirals in the two chains 
are opposite. Any finite J 12 will thus frustrate this spiral and the spins will prefer to point 
in the x direction so that J12 is unfrustrated. The classical ground state is then 

s^^iaoX-i)*, 

5 2 , i = 5 , (l,0,0)(-l)* +1 . (5.8) 



The linearized spin-wave spectrum about ( |5.8[ ) consists of two twofold degenerate 
branches with dispersions 

u^\q) = 2JS ^1 + 2Ji 2 / J — cos 2 (f) cos 2 q — sin 2 sin 2 q ± 2 cos cos q\j {J 12/ J) 2 + sin 2 sin 2 q 



In Fig ||a these two spin-wave dispersions are plotted for small q, J i2 = and = 0.02. As 
expected, in this case the energy spectrum is the same as for two isotropic spin chains, the 
only difference being that the zeros of the spectra have been shifted from q = to q = ±0 
(this shift occurs because £7(±0) does not commuting with the translation operator). Figure 
|5|b then shows u)(~'(q) for J12 = 0.01J. For these parameters the spectrum has evolved into 
a low lying 'acoustic' branch and a high energy 'optic' branch (not shown in the figure) in 
which the spins in different chains precess, respectively, in and out-of relative AFM phase. 
Note that the q = acoustic mode has acquired a gap because of the Ising-like isotropy 
induced by J 12 . 

The corresponding problem for two coupled planes is somewhat more complex. The 
effective Hamiltonian for a single CUO2-Y-CUO2 bilayer is again of the form 

H (2D) = H (2D) + H (2D) + H (2D ) (5JQ) 

where H[ 2D ^ describes the 'upper' layer in which the in-plane oxygen ions are buckled down- 
wards 

-^1 ^ = J y^\^i.i^\.i+-k + cos 0(^i^i,i+x + Si,iSt,i+x) + s i n 0(^M^i,i+x ~~ ^i,j^i,i+x) 

i 

+S'i,i^M+y + cos <fi(Si ti Sf ji+ 2 + S^jS^j+y) — sin0(5'f )i 5'i ii+y — S^Sf^)}, 

(5.11) 

describes the 'lower' layer in which the in-plane oxygen ions are buckled upwards 

H2 ^ — J^2{^2,i^2,i+x + COS + ^2,i^2,i+x) ~ S ^ n ^{^2,i^2,i+x ~ ^2,i^2,i+x) 

i 

+-S , 2 c >i 5'f j j + ^ + cos0(S , fjSf )i+y + 5 , 2 )i S , | j j +y ) + sin0(S'| j jS'| )i+y — 5 , | i j5'| )i+y )|, 



and we again assume that the layers are coupled by 

Hil D) = Ji2 y ZSi, i -S2, i . (5.13) 



The Hamiltonians ( |5.11|) and (|5.12 ) are obtained as before by taking the large U/t limit 



of ( |3.7| ) with Xij vectors given by ( |3.18| ) and ( |3.20| ). Here <p = \/t which we estimate to be 



0.02. A similar Hamiltonian (without the cos0 terms) was derived by Coffey, Rice and 
Zhang [] l| for a single Cu-0 layer in YBa2Cu306+ x - These authors noted that the classical 



ground state of a single layer is not a (possibly canted) Neel state. Instead, because 4>ij does 
not alternate in sign as in La2- x Nd x Cu04 there is no spin canting and the spins want to form 



a spiral [|TT],|T5|. It is interesting to note that the problem of finding the classical ground 
state of (|5.10|) is nontrivial because bonds which point in the x-direction favor spiraling 
about the |/-axis, while bonds which point in the y-direction favor spiraling about the x-axis 
and there is no classical ground state which spirals in this way. Fortunately, any difficulties 
associated with finding the classical ground state of ( |5.11| ) are probably irrelevant; just as 
for the coupled chains discussed above the spiral tendencies in the two layers are opposed to 
one another and so interlayer exchange, if sufficiently strong, will lock in a commensurate 
AFM state. 

Consider the set of classical Neel state parameterized as 

Si,i = >S(cos7 cos 77, sin 7 cos 77, sin?7)(— \^ Xi+y ^ 

$2,1 — 5'(cos7 cos?7, sin7 cos?7, sin?7)(— l)^ i+?/i+1 ' ) . (5-14) 



Treating these as variational states for ( |5.10| ) yields an energy per site of 



E[<f>,rj\ = -JS 2 (cos(f)(sm 2 r] + l) + cos 2 ri) - J l2 S 2 (5.15) 

which is minimized when 77 = 0. Thus the lowest energy Neel states are those in which the 
spins lie in the xy plane. The origin of this 'easy plane' is rather subtle. It arises from 
the cos (pij easy-axis terms in ( |5.10| ). These terms favor spins parallel to the x direction on 



bonds which point in the y direction, and visa versa, so that when a classical Neel state is 
rotated in the plane the energy gain on one type of bond increases while that on the other 
type decreases in such a way that the total energy remains constant. 

The classical spin- wave spectrum about (|5.14 ) when rj = can again be calculated using 



standard methods. The resulting spectrum has four branches with dispersions given by 
#(?) =2JS( { A(0) + 2J l2 /J)A(0) - A(q)B(q) - C(q} 2 



± ((A(0)B(q) - (A(0) + 2J 12 /J)A(q)) 2 + AC(q) 2 A(q)B(q)] 

(5.16) 

^(q) = 2JS^( A(0) + 2J 12 /J)A(0) - A(q)B(q) - C(q) 2 

K 1 

± (((A(0) + 2J 12 /J)B(q) - A(0)A(q)) 2 + AC (q) 2 A(q) B (q) 



where 



A(q) = (cos 2 7 + sin 2 7 cos </>) cos q x + (sin 2 7 + cos 2 7 cos 0) cos q y , 
B(q) = cos0(cosg :c + cosq y ), 

C(q) = sin 0(sin 7 sin q x + cos 7 sin q y ) . (5.18) 



It is now apparent that the Neel state ( |5.14| ) with r\ = is not stable for small J\2j 'J r . In 



the limit q x ,q y ,(f) <C 1, Ju/J < 2 , and 7 = the dispersion of the lowest lying branch is 
well approximated by 

W H(g) ~ 2J5I I 7" 1 C 2^ i . (5.19) 




When J12/ J < 3<f) 2 / 4 the frequencies of this branch become complex for a region in g-space 
near q = 0. These complex frequencies signal the appearance of unstable modes which grow 
exponentially with time. 

In Fig. ^| the spin- wave dispersions given by ( |5.17|) and (|5.18|) are shown for 7 = 0, q y = 
and q x small when = 0.02 and J-&/J = 0.0, 2 /2, 30 2 /4 and 0.01. The hatched regions 
denote q values for which the spin wave frequencies become complex indicating that the Neel 
state is unstable. As J i2 is increased the spin-wave spectrum evolves in the following way. 
First, for J 12 = 0.0 there are two twofold degenerate branches, the lower two of which become 
complex at small q x . Then, as J12/J is gradually increased the degenerate branches split 
and the region of instability in g-space shrinks so that when Ju/J — <p 2 /2 there is only one 
unstable branch (u(~') and when Ju/J — 3</> 2 /4 the complex frequencies disappear entirely 
from the spectrum and the Neel state becomes locally, and almost certainly globally, stable. 

For the physically relevant case J ^> J 12 <f) 2 J the dispersion of the four spin- wave 
branches are approximately given by 



U3 



(") 



(q) ~ 2V2 JS\q\ (5.20) 
J b -\q) ~ 2V2JS^<J) 2 + \q\ 2 (5.21) 
J a + \q) ~ ul + \q) ~ 2JSy/BJ 12 /J+\qi 2 . (5-22) 

Such a spin- wave spectrum is precisely what is expected for a bilayer made up of two easy- 
plane antiferromagnets [^(J. In this limit uj^~\q) and uo[ \q) have evolved into low energy 
modes corresponding to gapless in-plane and gapped out-of-plane spin-waves, respectively. 
These modes are 'acoustic' modes in the sense defined above: the spins in the two layers 
precess in relative AFM phase with one another. Figure |6|d shows these low lying modes for 
J 12 / J ~ 0.01. The out-of-plane mode has a gap of 2\^2JS(j) a result which is consistent with 
the experimental observation of a gapped out-of-plane zone-center mode 0. The remaining 
branches u^\q) and u^ + \q) have evolved into the high-energy 'optic' modes in which spins 
in the two layers precess out of AFM phase and have a gap of 4v / 25'v / J12J ■ These modes have 
not been seen experimentally for energies up to ~ 50 meV indicating that J i2 > 0.01J. 
The interlayer coupling is therefore large enough to lock in the Neel state according to the 
above scenario. 



VI. CONCLUSIONS 



In this paper the anisotropic corrections to superexchange arising from SO coupling in 
the distorted Cu-0 layers in La^^Nd^CuC^ (tilting distortion) and YBa 2 Cu 3 6+:1 . (buckling 
distortion) have been calculated using Moriya's method Special care has been taken to 
include the higher-order symmetric anisotropy terms whose importance has recently been 
emphasized by SEA ||. 

In the La2- a; Nd a ;Cu04 system in the presence of a tilting distortion it was shown that 
regardless of the tilt axis it is always possible to find a local coordinate system in spin space 
in which the SO induced spin precession of moving electrons has the same form. This form is 
precisely that which was studied in J?| and is described by Hamiltonian In this special 

spin space coordinate system the classical ground state of the effective spin Hamiltonian has 
no weak ferromagnetic moment but when the system is viewed in the physical coordinate 
system a weak ferromagnetic moment appears. The ratio of the octahedral tilt angle to the 
weak ferromagnetic moment depends on the tilt axis according to (|4.6| ) and so it decreases 
as the tilt axis moves from (110) in the LTO phase to (100) in a hypothetical insulating 
AFM LTT material. The experimental observation of the reduction of this ratio in the Pccn 
phase of La2- a; Nd a ;Cu04 would provide a test of the theory discussed here and in ||. One 
unresolved problem with the theory is that it predicts that regardless of the tilt axis there 
should be only one spin-wave gap - a prediction which is at odds with the experimental 
observation of two different gaps for the in-plane and out-of-plane zone-center modes. It is 
likely that this discrepancy is due to the additional anisotropy caused by dipolar interactions 
between the spins. However, a full treatment of the problem including these interactions 
has not yet been carried out. 

When the same method for calculating the SO corrections to superexchange in 
La 2 _ x Nd x Cu04 is applied to the buckled CUO2-Y-CUO2 bilayers in YBa2Cus06+x the re- 
sulting Hamiltonian is quite different from what has been used previously to model this 
system. The Hamiltonian obtained here describes a system in which the spins in each indi- 
vidual Cu-0 layer tend to form a spiral pattern. However, one consequence of the inversion 
symmetry of the unit cell of YBa2Cu306+ x is that the senses of the spirals favored by the 
two Cu-0 planes in a given bilayer are opposed to one another. When a weak interlayer 
coupling is included the spiraling become frustrated and a commensurate AFM state is sta- 
bilized. For physical parameters the low-energy spin-wave spectrum agrees with what is 
seen experimentally: there is an in-plane gapless mode and an out-of-plane gapped mode. 

We conclude on a speculative note. In |7j it was argued that the frustrated DM terms in 
La2Cu04 can give rise to new physics away from half filling. In particular, doped holes can 
gain energy from their SO induced spin precession as they move through a commensurate 
spin background. This energy gain is similar to that which occurs when holes move through a 
spiral state [0]. We have shown here that frustrated DM terms are present in YBa 2 Cu30 6+ . r , 
and the size of these terms is significantly larger than in La2Cu04 (by roughly a factor of 
five). Because of this, holes doped into a Cu02-Y-Cu0 2 bilayer can also gain energy through 
their spin precession. In this case the maximum energy gain occurs when the spins form a 
commensurate Neel background and are lined up parallel to the z direction. When moving 
through such a spin background doped holes see, effectively, a double spiral ^7|. We note 



that a strong tendency for commensurate spin fluctuations in superconducting YBa 2 Cu306+ x 



has been observed in both NMR and neutron scattering P,pT| experiments while in doped 
La2Cu04 spin fluctuations appear to be incommensurate |29|j30[| . It is tempting to speculate 
that SO effects such as those discussed here and in [[7| may play some role in determining 
the different spin dynamics of these two systems. 
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In this appendix the exact unitary transformation which maps Hamiltonian (|3.7|) into 
( |1.1[ ) is constructed for arbitrary \ an d 0- This transformation is most easily represented as 
a sequence of two local rotations in spin space. The first rotation eliminates the component 
of the unfrustrated precession perpendicular to the frustrated axis and is generated by 



(1) 



U (1) = exp jz (-1)^+^ tan- 1 ^ 2 S ™ 2X )(cosx, - sin X , 0) • S, 

When is applied to Hamiltonian fl3.7| ) the result is 

C/(i) OT (i)t =H' = J2 \ c l (v / * 2 + («2Sin2 X ) 2 ^ + ■ a aP ) c jP + h.c 

m 1 v J 

a/3 

+uJ2 n n n n ( 2 ) 



where the transformed A' vectors are 

K,i+* = (-l) Xi+Vi (ai - "2 cos 2x) (cos x, sin x, 0), 

K A+ y = -(-l) !Bi4lw (ai + a 2 cos2x)(cosx,sinx,0). (3) 

The precession axes on all the links in the lattice are now parallel to one another. A further 
rotation generated by 

= exp (^(-l^ W 1 ( a * COs2 X ) ( C os X , sin X , 0) • s\ . (4) 

k 2 + a!sin 2 2x/ J 



then ensures that the precession has the same magnitude on all the links. Upon applying 
this transformation the Hamiltonian becomes 

UWH'UW = £ { c f Q fe Safj + £». . B aP ) c jP + h.c.} + UJ2 n^n n (5) 

(i.n\ i 



where £^4* = Jt 2 + a 2 , + aia2sin2x, £i,j+ y = Jt 2 + a^ — aia2sin2x; and 



A-, + x = (-l)* i+% ai (cos X , sin X ,0) (6) 

If the second order in 9 corrections to the direct hopping integral are ignored (£y ~ t) and 
a global rotation in spin space is performed so that the z-axis is parallel to (cosx, sinx, 0) 
the result is Hamiltonian flO]) with M+X = (-l^+w) taa~ 1 (a 1 /t) ~ (-1) (Xi+k) 0.05# and 
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FIGURES 



FIG. 1. Copper xz orbitals and oxygen a orbital in a Cu-O-Cu bond parallel to the x di- 
rection in the presence of (a) a tilting distortion (La2_ a; Nd a ;Cu04) and (b) a buckling distortion 
(YBa2Cu306+x)- In (a) the the xz orbitals are tilted with respect to the copper-oxide plane while 
in (b) they are not. As a result, the ratio of the Cu-0 bond angle to t(xz, a) (the hopping am- 
plitude between the xz and a orbitals) is larger by a factor of ~ 5 in YBa2Cu306+a; than in 
La2- a; Nd a ;Cu04. 

FIG. 2. Unfrustrated (a) and frustrated (b) spin precession about a single plaquette in the 
presence of a tilting distortion. The precession shown in (a) is unfrustrated because the electron 
spin returns to its original value upon hopping around any closed loop. The precession shown 
in (b) is frustrated because an electron which hops around a single plaquette acquires a finite 
spin precession. While it is always possible to eliminate the unfrustrated precession by a unitary 
transformation the frustrated precession cannot be so removed and thus is responsible for any 
physical anisotropy in spin space due to spin-orbit coupling. 

FIG. 3. Two neighboring copper-oxide clusters in the Pccn phase of La2- ;c Nd; ! ;Cu04. The 
heavy line is the (cos x> sinx, 0) axis about which the CuOq octahedra are tilted through an angle 
9. The light line is the (cos x, — sinx, 0) axis. The spin precession of a single electron which hops 
in the presence of this tilting distortion is a combination of frustrated precession about the light 
axis and unfrustrated precession about the heavy axis. The classical ground state of the effective 
spin Hamiltonian describing this system at half filling is also shown. The spins are nearly lined up 
along the frustrated axis except for a slight cant out of the plane by an angle which is proportional 
to sin(2x). 

FIG. 4. Two buckled CuO 'chains' representing the toy model described by ( |5.1| ). The spin 
configuration shown is one possible classical ground state in the absence of interchain coupling. 
The DM interactions induced by the buckling cause the spins to form a spiral, but because the 
buckling in the bottom chain is opposite to that of the top chain the spirals in the two chains have 
opposite senses. 

FIG. 5. Spin-wave spectrum for two buckled chains when (a) J\2 = and (b) J%% = 0.01J. In 
(a) the spin wave spectrum is identical to that of two isotropic Heisenberg models except that q 
has been shifted by ±</>. In (b) there is a finite interlayer coupling and the system has acquired an 
Ising-like character and a gap. 



FIG. 6. Spin-wave spectrum about a Neel state with spins parallel to the x direction for a 
Cu0 2 -Y-Cu0 2 bilayer in YBa 2 Cu 3 6+;E for J 12 /J equal to (a) 0; (b) 0.25</>; (c) 0.75^; and (d) 0.01 
where (j) ~ 0.02. The hatched regions represent q values where one or more of the frequencies are 
complex indicating an instability. There are four modes which degenerate into two when Ji 2 / J = 
(a). As Jn/J increases these modes split. In (b) one of the two unstable modes has become stable. 
In (c) both have become stable and in (d) the two higher lying 'optic' modes are not shown and the 
two lowest modes correspond to a gapless in-plane and a gapped out-of-plane mode in agreement 
with experiment. 



